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Outage Probability Analysis of Dual Hop 
Relay Networks in Presence of Interference 

Bappi Barua, Mehran Abolhasan, Daniel Franklin and Farzad Safaei 



Abstract 



Cooperative relaying improves the performance of wireless networks by forming a network of 
multiple independent virtual sources transmitting the same information as the source node. However, 



interference induced in the network reduces the performance of cooperative communications. In this 
work the statistical properties, the cumulative distribution function (CDF) and the probability density 
function (PDF) for a basic dual hop cooperative relay network with an arbitrary number of interferers 
over Rayleigh fading channels are derived. Two system models are considered: in the first system model, 
the interferers are only at the relay node; and in the second system model, interferers are both at the 
relay and the destination. This work is further extended to Nakagami-rn faded interfering channels. 
Simulation results are presented on outage probability performance to verify the theoretical analysis. 

Index Terms 

Signal to interference plus noise ratio, interference to noise ratio, amplify and forward, probability 
density function, cumulative distribution function. 
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I. Introduction 



Wireless channels often suffer severe performance degradation due to multipath fading and 
interference. The increase of spectral reuse in wireless systems exposes the network to ever- 
large number of interfering nodes. Ad-hoc and multiuser networks often face interference from 
other user nodes operating in multiple parallel channels. Interference can also appear in the 
network as cochannel interference. Cooperative relay communications, originally introduced by 
van der Meulen [1] has been extensively studied in CD-[[51 without the presence of interference. 
In cooperative relay networks, cooperating nodes help the information source node by jointly 
forming a virtual antenna array and thus providing higher diversity gain and information reception 
reliability. In addition, cooperating relays can offer a very flexible extension of network coverage 
without substantial development of infrastructure. However to date, research efforts on the 
interference channels have mainly concentrated on information theoretic analysis lloll- lfTTl and 
using assisted relay networks in lfT2l - ll26l . Statistical modeling for single hop multiple antenna 
systems with active interferers has been performed for maximal ratio combining (MRC) receivers 
in ll27l - ll3Tl and using optimum combining receivers in Il32l - ll35l . 

Interference in cooperative relay networks has been extensively studied in previous publications 
[fl2l - [|22l . Authors of llT2l and lfT3l consider a scenario where only the relay node experiences 
interference in interference limited networks and the receiver node remains interference free. The 
relays here can estimate the instantaneous CSI of interfering channels to scale the gain. A fixed 
gain AF relay with interference limited destination is considered in Ifl4l . However, the assumption 
that the AF relay gain parameter includes the instantaneous or average channel information of 
interfering channels (as assumed in lfl2l - lfT5l . [TTTll . [[191 , [1201 ) requires additional computational 
capability at the relaying node, and in certain cases where the interfering signals are not known 
to the relay a priori, the technique can not be applied. Furthermore, in many previous results 
authors consider interference limited cooperative networks lfT71 - [|22l . Outage performance of 
a dual hop network has been studied using a fixed gain relay in ifTTl and hypothetical gain 
AF relay in [fT8l with an arbitrary number of interferers. The system is assumed to be an 
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interference dominated network where noise power is negligible compared to interfering signal 
power. Closed form PDF and CDF of signal to interference ration (SIR) are derived. Avoiding the 
noise in receiver systems usually provides an opportunity to express some of the integrals used 
to obtain the CDF and the PDF in very standard formats. Outage probability using DF relays 
in Nakagami-m fading channels is studied in ll22ll . The authors consider that the destination 
faces a negligible amount of interference, but non-negligible noise. The PDF and CDF of 
corresponding receiver SNR is derived. In ||25ll , the authors have studied transmission techniques 
designed to manage interference using an information-theoretic approach. Bit error rate (BER) 
performance has been presented in ll36ll for a multi-user network in which multiple sources 
communicate with one receiver node via a common relay in a dual hop network. A different 
approach has also been proposed to cancel the interference at the relay node. In [37], the authors 
have investigated the effect of cooperation in an interference-limited system utilizing a decode- 
and-forward (DF) relaying protocol over Rayleigh fading channels. It divides the cooperating 
nodes into different cooperating regions and quantifies the relation between cooperative region 
radius and the interference level. It also analyzes the network sum rate optimization problem for 
a cooperative region radius. 

In this work, we derive the closed form expressions of statistical properties, the CDF and 
the PDF of signal to noise plus interference ratio (SINR) of an interference relay network that 
is operating under the influence of arbitrary number of interferes and noise. An amplify-and- 
forward (AF) relaying protocol is adopted due to its simplicity and ease of deployment at relay 
nodes. Lastly numerical results on the outage probability performance are presented to verify 
the theoretical analysis via Monte-Carlo simulations. 

II. System and Channel Model 

Consider a single source-destination pair communicating via a single antenna relay without any 
direct link. We will denote source-relay and relay-destination links as S-R and R-D respectively. 
A half duplex AF protocol has been considered over independent and non identically distributed 
(i.n.i.d.) Rayleigh fading channels. In this paper, two different system models are investigated: 



system model 1 (SM 1), in which the interferers are only at the relay node; and system model 
2 (SM 2), where interferers affect both the relay and destination nodes. In both models, all the 
interfering channels are i.n.i.d. Rayleigh faded. The destination is assumed to have full channel 
state information (CSI) of the two main channels, S-R and R-D, while the relay has full CSI of 
the S-R channel only. The source and relay have no CSI of forwarding transmitting channels. 
None of the nodes, source (S), relay (R) and destination (D) possess information about the 
interfering channels. 

The instantaneous and average signal to noise ratio (SNR) of 1st and 2nd hops are denoted as 
7i — P ^2 and Aj = ^r- respectively, where i G {1, 2}, P is the corresponding source and relay 

i i 

power; hi and £1^ are the instantaneous and average channel gain of the ith hop respectively^ and 
of is the one sided additive white Gaussian noise (AWGN) power at relay or destination node, 
i.e. % G {R, D}. We assume the total power of the network (source and relay) is constrained to 
P tot . This total power is split between the source and the relay by a power sharing coefficient 
C G (0, 1] such that the source and the relay powers are given by, P s = (1— C)-Ptot and P R = CAot 
respectively. Thus if equal power sharing protocol is adopted, ( will be 0.5. 

Let there be a total of L interferers in the system, and define an interferer set 3, the set of all 

interfering source nodes. For example, any interferer // G 3, where I E C, C = {1,2,3 ,L}. 

The elements of interfering channel row vectors hf « G C L represent the corresponding interfer- 
ence channels from the source element nodes of the interfering signal source vectors x/ ti G C L , 
i G {R, D}u Furthermore, we assume the interferer Ii has individual transmit power Pj Thus 
the instantaneous and the average interference power to noise ratio (INR) for any interferer // 

I i2 

is — / an d Aj^ = p '^ 1 ' 1 respectively, hi i and af are the fading channel gain from 
the interfering source Ii to the node i and the noise power at node i respectively, / G C and 
% G {R, D}. Qjj is the average interfering channel gain, Qjj = E [|/iz,i| 2 ] . Throughout this paper 
the term INR will be used to indicate individual interferer signal power to noise ratio unless 



'The average channel gain, f2,: is in fact the statistical average of the squared instantaneous channel gain hi, i.e. fi; = E [|/ii| 2 ] . 
2 C L denotes a L-dimensional complex vector. 
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A. System Model 1 : Interference at Relay 

In the network model as shown in Fig. [H interference occurs only at the relay node. The 
received signal at the relay node R is 

y R = hix s + h/^xf^ + n x (1) 

where, hfi G C L and x/ i G C L are the row vectors for the interference channels and the 
corresponding interfering signal sources respectively. The notation (.) T represents the transpose 
operation on a vector. The signal at the destination node is 

y D = Gh 2 h 1 x s + Gh 2 h Itl xJ A + Gh 2 n 1 + n 2 (2) 

where G is the AF relay gain, rii ~ CM (0, af) and n 2 ~ CM (0, a 2 ) are additive white Gaussian 
noise at the relay and destination respectively!] At the receiver node, the signal to interference 
plus noise ratio (SINR) will be 

G^h^h^Ps 

7SINRl G^lXxS/X, + G*\h 2 \ 2 a! + al {) 

where the diagonal matrix £/ = E |x| jX/,1 j is composed of interference signal powers. With 
hypothetical AF relay gain G 2 = p ^a > tne end-to-end SINR is 

7172 fA\ 

7siNRi = ; ; (4) 

7i + 72 + 727/ 

where 71 and 7 2 are the instantaneous SNRs for S-R and R-D links respectively. 7/ is the total 
interference to noise ratio at the relay node, which is literally the sum of all individual INRs. 

Note that due to the Rayleigh fading assumption, the first and second hop SNR 7x and 72 are 
exponentially distributed with mean parameter Ai and A2 respectively. Later, it will be shown that 
under an equal power sharing protocol at the source and relay nodes, the system performance 

3 CA/"(A,cr 2 ) denotes a circularly symmetric complex Gaussian random variable with mean A and variance a 2 . 
4 A CSI assisted relay gain G 2 = ^ +(y >± proposed in (2) render the end-to-end SINR as 7sinr = 71+7 ^_?^ 7f +1 ■ Eq- © 
proposes a tight upper bound on the CSI assisted SINR 7sinr. 
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metric of system model 1 does not vary if the interferers are switched from the relay to the 
destination node. 



B. System Model 2: Interference at Relay and Destination 

In this network model, interference occurs both at the relay and at the destination nodes. 
We let there are L x interferers at the relay node and L 2 interferers at the destination. All the 
interfering channels are independent and non-identically distributed. In this case, the received 
signal at the destination is 

y D = Gh 2 h 1 x s + G , /i 2 h/ i ixJ )1 + h/ j2 x^ 2 + Gh 2 n\ + n 2 (5) 

where xj ^ € C Ll and x /2 £ C L2 are the vectors with interference sources for the relay and 
receiver nodes respectively, and h/i e C Ll and h/ 2 € C Li are the corresponding fading channels 
from interferers to the relay and receiver nodes respectively. Again, as for SM 1, G is the AF 
relay gain, n x ~ CM (0, of) and n 2 ~ CM (0, of) are AWGN at the relay and the destination 
respectively. Thus the SINR with arbitrary relay gain is given by 

7SINRa G^lXiS/Xi + h i^iM,2 + G 2 \h 2 \ 2 aj + a\ 

S^i = E |x| ^x/^j and S/ j2 = £ |x} 2 x/ )2 | are diagonal matrices of the transmission powers 
of interfering signals at the relay node and the destination respectively. Applying the similar 
hypothetical relay gain as used in SM 1 at the AF relay, the receiver SINR is given by 

7i72 (n . 

7SINR 2 = ; ; ; (7) 

7i + 72 + 7i7/,2 + 727/,i 
where 7 /;1 and 7/ 2 are the total interference to noise ratio at the relay node and the destination 
respectively. 



III. Statistical Analysis 
In the following sections statistical theorems related to the network of interest are developed. 
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Definition 1: Let random variables X, Y and U be statistically independent, where X 
£ {j^j an d Y ~ £ (j^J J3 Random variable U — J2 Ui + 1, where Ui ~ £ (x - ) are i-n.i.d. expo 



C(A)r»(A) 



nentially distributed random variables and U has the PDF, fy {u) = Yl Yl Xl,J ^^ u j 1 - > — e 

i=i j=i r w A «<i> 

u > 1 OS. A r.v. W is defined as 



n 



W = XY (8) 

Theorem 1 (CDF and PDF): The cumulative distribution function (CDF) and probability den- 
sity function (PDF) of random variable W are given respectively as 



C(A) Ti(A) 

F w («,) = 1 - £ jfy (A)e" 

i=l j = l 



A, 



oo k+1 

*+EE 



r (j) (X x Xy) + k\(k + 1)! VA + w\ u 



X x + wX u (i) 

2T (j +n)Un 



k=0 n=0 



k + 1 

n 



w 



\/X X Xy 



+ G 



1,3 
.3.2 



/ X X X U ^ 


l-J-n,l,l\ 


\X X + w\ u (i) 


1,0 / 



(9) 



C(A)7i(A) 

fw h = E E x « ( A ) e 
i=i j=i 

OO fc + l ^ 

-EE 



An 



As + wX u (i) 



k=0 n=0 



w 



2fc+l 



T(j)(X x X y ) k+1 k\(k + l)\ \X x + wX u{l) 



X x X u (i) 



_A X Ay A x + wX u (i) 

2T (j + n) + Zi (w, A) <( 2A; + 2 



tw (A* + Ay) (j + n) wA„(i) I wX u{i) ,, : . 

+ ~ — : : ^4,3 



X x X y 



X x + wX u {i) X x + wX u (i) 



X x X u (i) 
X x + wX u (i) 



1-j-n, 1,1,0 
1,1,0 



(10) 



5 £ (x - ) denotes an exponential distribution with a hazard rate j-. 

6 where A = diag (A U1 , \ U2 , A„ L ). £ (A) is the number of distinct diagonal elements of A and X u n\ > A u ( 2 > > > 

Au(c(A)) are the distinct diagonal elements in decreasing order, n (A) is the multiplicity of \ u <i) and Xi t j (A) is the (i, j')th 
characteristic coefficient of A 1381 . 



where ^ = \ V (* + 1) + \ 1> (k + 2), and J x (u;,A) ^ G$ ^5^- 



1 - J - 71, 1, 1 

1,0 



+ 



Up 



6, 



and ijj (x) are the Meijer-G function and the 



Euler psi function defined as ll40l eq. 8.2.1.1] and PT1 eq. 8.360.1] respectively. 

Note that, exchanging X and Y in Definition 1, represents the switching of the interference 
from relay to destination node in ©. Hence, the CDF and PDF equations of SMI are also 
applicable to the general case when either the relay or destination is subject to interference. In 
addition, the performance metric will be the same if an equal power sharing is adopted between 
the source and relay nodes, that is, when ( = 0.5. 

Proof: See Appendix [A] ■ 

Definition 2: Let random variables X, Y, U and V are statistically independent where X ~ 

£ and Y ~ £ (■}-). Random variables U = £ U t + 1 and V = E v i + 1. u i ^ d v i are 

^ ^ y ' 1=1 1=1 

i.n.i.d. exponentially distributed random variables with PDFs ll39l 



C(Ai)T t (A0 (u _ 1} 

i=i j=i 1 UJ A «{i) 



C(A 2 )r,(A 2 ) 

^ ^ r (g) A* (p) 



respectively. We define a r.v. Z such that 



Z = Xy (13) 
XU + YV 



Theorem 2 (CDF and PDF): The cumulative distribution function (CDF) and probability den- 
sity function (PDF) of random variable Z are given respectively as 
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C(A 1 )r 8 (A 1 )C(A 2 )r p (A 2 ) 

^W = 1 "EEEE X m (Ai) X M (A 2 ) e 

i=l j = l p=l q=l 

oo fc+1 fc+1 



X, 



Ax 



x 



X x + z\ v (p) 
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x G 
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(14) 
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Proof: See Appendix [B 
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Proposition 1: Consider a system with L interferers in both relay and destination nodes where 
the interfering channels are i.i.d. Rayleigh faded. Adopting an equal power allocation protocol 
for source and relay node the cumulative distribution function of 7sinr 2 in © can be written as 



F (V) = 1 - e" 27/A ( - 



2L 



T(L)' 



-2 7 /a 



00 fe+i fe+i 

k=0 n=0 r=0 



k + l\ (k + 1 
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1,0 



(16) 



Proof: See Appendix | 



IV. Interference at Relay: I.I.D. NAKAGAMi-m Interferers 

Proposition 2: Suppose the interfering channels are i.i.d. Nakagami-m distributed while the 
main channels (S-R and R-D) are Rayleigh faded. In this case the CDF of 7sinr 2 can be written 

as 



7SINR 2 



( 7 ) = l-e" 
00 fc+i fc+i 

k=0 

Aiaj,2 
Ai + 7 a/, 2 



7 
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(17) 
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where ai t i = Xi^/rrii, i E {1, 2}. A/,1 and A/^ are the average INR at the relay and the destination 
respectively; similarly, mi and rri2 are i.i.d. Nakagami-m parameter at the relay and destination 
respectively. The average SNR of the 1st and 2nd hop are Ai and A 2 respectively. 

Proof: See Appendix iDl ■ 

V. Outage Probability 

In this section the derived analytical results are used to the investigate wireless network 
outage probability performance. We define outage probability as the probability that the in- 
stantaneous receiver SINR falls below a predefined threshold value of SINR 7th. We consider 
7th = P (2 MR — l), where p varies from 1 to 6.4 depending on the degree of coding, M be the 
number of hops and R the data rate in bits/s/Hz 02). The CDF equations ©, (O and (fTTT) 
may be used to evaluate the outage probabilities in SM 1 and SM 2 with an arbitrary number 
of interferers and interfering powers. 

VI. Numerical Analysis 

In this section, numerical results on the derived equations are presented and compared with 
Monte-Carlo simulations. Due to lack of transmitter CSI we assume the source and the relay 
evenly share the total system power P tot . Furthermore, for fair comparison with relevant studies 
the average channel gain of all hops in this network are set to unity and the noise variance at 
the relay and destination nodes are set to a 2 . We consider 7 th = 3 which achieve a data rate 
R — 1 in a dual hop network unless otherwise specified. The maximum value of k = 100 for 
sum terms in © and (fl4l) is used, which is sufficient to match the analytical result with the 
simulation data. The number of interferers at the relay and destination are assumed to be same 
(Li = L 2 ) in SM 2 networks, and all the interferers have equal transmit power. Throughout this 
section we will use total average transmit SNR, A to t — -Ptot/c 2 as a function in performance of 
SINR plots. 

Fig. [2] shows the outage probability (OP) as a function of total transmit SNR for a total of 
{4, 8, 16} interferers in SM 1 and SM 2 with INR fixed at 3 dB. All wireless channels (main 
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channels and the interfering channels) are assumed to be Rayleigh faded. As expected, outage 
probability increases with an increase in the number of interferers. However, when SM 1 has 
the same number of interferers as the total number of interferers for SM 2, SMI shows slightly 
better performance compared to SM2 at a moderately low SINR region only. To study the OP 
performance with different interfering powers Fig. [3] shows outage probability as a function of 
A to t for a system with 4 interferers at the relay and 4 interferers at the destination. The INR 
values {3, 6, 9, 12} dB are considered. Interestingly, the figure suggests that every 3 dB increase 
in average INR values requires a subsequent increase of 3 dB total SNR on average to ensure 
an OP of around 1CT 3 . In contrast, Fig. @] shows OP in a fixed 3 dB INR system in SM 2 with 
a total of {4, 8, 12, 16} number of interferers at the relay and destination respectively. To see 
how the interference system behaves if one of the nodes experiences most of the interference, 
Fig. H] shows OP with a large L 2 and L\ = 4. For very high interference levels e.g. L 2 = 10000 
interferers at the destination node, a A tot = 77 dB is required to achieve an outage probability 
of 10~ 3 compared to 50dB when L 2 = 4. 

Fig. [5] shows outage probability with number of interferers L in relay and destinations where 
Li = L 2 = L and each interferer has an INR of 3dB. The figure shows how the outage 
probability increases exponentially as the number of interferers increases. The plot includes 
graphs for transmit SNRs of 20, 30, 40, 50 and 60 dB. A very interesting phenomenon observed 
is that at 60dB of total transmit SNR the system can support up to 40 interferers at the relay 
and destination, ensuring an outage probability of 10~ 3 , however, with a reduction of only lOdB 
in total transmit SNR the performance reduces drastically and it can support only 4 interferers 
at the relay and destination for the same outage probability. 

Fig. [6] shows the outage probability as a function of interference power when the relay 
and destination are both subject to 4 interferers. The figure compares the outage probability 
of an interference limited system (dashed line) studied in [fT8l to a system where both noise 
and interference are present. To express total transmit power and interference power in dB we 
normalize the terms considering noise power equal to 1 unit. Total transmission power of 20, 
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30 and 40 dB are considered in the figure and the power of each interferer is rated from -30dB 
to +30dB. Figure shows that the outage probability the noise plus interference network will be 
similar to an interference limited network when the transmission power of each interferer is more 
than 10 dB. However, performance will vary if the number of interferers in the network changes. 
The figure also illustrates how noise dominates outage probability metric when interference power 
is low. 

The PDF of 7sinr 2 is presented in Fig. [7] for two different A tot , 20 dB and 25 dB, when the 
individual interferer INRs are {3, 6, 9} dB. The figure implies that a lower INR per interferer 
increases the probability of higher output SINR at the receiver and similarly when the INR is 
constant, an increase in average total signal power results in a higher output SINR. In addition, 
Fig. [8] shows the CDF of 7sinr 2 as a function of 7 th in dB for two A tot groups, 20 dB and 30 
dB. Each group contains plots for 3 INRs 3, 6 and 9 dB. It reveals that higher average total 
SNR and lower INR per interferer decreases the probability of outage. 

Finally, Fig. [9] presents a plot over Nakagami-m faded interfering channels. The outage 
probability is plotted as a function of A to t for different number of interferers and different 
Nakagami channel parameter m in SM 2. Surprisingly, varying the Nakagami m parameter does 
not result in any significant effects on outage performance if the number of interferers and INR 
remain constant in the network. 

VII. Conclusion 

In this work, we derive the cumulative distribution function and the probability density function 
of a dual hop interference relay network with arbitrary number of interferers, where the main 
channels and the interfering channels are i.n.i.d. Rayleigh faded. Later the analysis is extended 
for i.i.d. Nakagami-m faded interfering channels. This analysis has been performed for two 
system models; in SM 1 interference is only at the relay node while in SM 2, interference 
is at both the relay and destination nodes. The derived CDF and PDF expressions for SM 1 
can also be used for analyzing an interference network where interference occurs only at the 
destination node. Numerical results for outage probability performance for different network 
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parameter configurations are presented. The Monte-Carlo simulations show an exact match with 
the analytical expressions. Results indicate that in a total INR constrained interference network, 
the number of interferers do not affect the performance of the system in a large scale. Similarly, 
if the interfering channels are Nakagami-m distributed, changing in Nakagami-parameter m does 
not affect the performance of the system if the number of interferers and the INR per interferer 
remain constant. 



Appendix 

A. Proof of Theorem [7J 

According to the PDF of X, Y and U as defined in Definition 1, the random variables X, Y 
and U are nonnegative, thus F w (w) = for w < 0. For w>0we have 

r XY 

= P {X (Y - w) < wYU} 

wyu 

poo poo r~_ — 

= 1 - e~ w/x « + / / I fx 0) h (V) fu (u) dxdydu (18) 
Ji i» Jo 

The two inner integrals in (fT8l can easily be solved by using BTl eq. 3.351.3] and Pffl 
eq. 3.471.9] and after some manipulation the CDF of W is obtained, 



C(A)r l (A) , . p -w/\ v roo 

Fw w = i - E E S /ax / Mu ~ 1 ) < " 1 ^ ( - 1,/i • , " 



e '™" x ' Ki ( 2w \[S) du (19) 



where K v (x) is the z/th order modified Bessel function of the second kind. Using pTl eq. 8.446] 
we expand the first order modified Bessel function of second kind, then HOl eq. 8.4.6.5], 
ll43l eq. 3.40.1.1] and Pffl eq. 3.351.3] are applied to solve the related integrals. After some 
manipulation we arrive at the desired result ©. The PDF of W follows directly from the 
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differentiation of the CDF of W in eq. © w.r.t. w by using EQl eq. 8.2.2.30], 

- _ T <r-ln m + 1 > n I r 



A. 

dx 



% g I X 





a p , -a \ 


^X 






l-a,b q J 



(20) 



B. Proof of Theorem \2\ 

Following a similar argument as used in the proof of Theorem CD the random variables X, Y, 
U and V are nonnegative, and thus F z (z) = for z < 0. For z > 0, 



F z (z) = P 



XY 



< z 



XU + YV 
¥{X(Y -zU) < zYV} 

fy (y) fu (u) dydu + 



oo ruz 



1 JO 



OO POO POO 



1 Jl Juz JO 



(21) 



fx (x) f Y (y) fv {v) fu (u) dxdydvdu 



A T 

— iuu 



V 

- T 

— J-uvyx 



The first part of this expression, integral I uy can be solved by using PTl eq. 3.351.3], resulting 
in 



(22) 



C(Ai) n (Ai) , . j 

In Ju^ys, the first two integrals for variables x and y can be solved using a similar approach, 
via EH eq. 3.351.3], EH eq. 3.471.9], [41, eq. 8.446], ffl, eq. 8.4.6.5] and EH eq. 3.40.1.1]. 



Substituting (|22l) in (T2TI) and after some manipulation we have the CDF of Z 



C(A 1 )r i (Ai)C(A 2 )r p (A 2 ) 



X M (AQX p , g (A 2 ) 

»=i j=i P =i q =\ r(j)r(g)A^A^ 
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EE 
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n 
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X x X v (p) 
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q+n 



^x^v(p) 
X x + Z\ v (p) 




u k+l (u - 1) 
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k=0 n=0 



n 



k\ (k + 1)! (X x X y ) k+1 \X x + zX 
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(u - e A «« e~ uz/x « \\nj -^—z -ip k \du 

X X Xy 



(23) 
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where ip k = \ ip (k + 1) + \ ip {k + 2). Now using [41, eq. 3.351.3], EH eq. 8.4.6.5] and gl 
eq. 3.40.1.1] it is possible to solve (l23l) . which then results in the desired eq. (fl"4l) . Eq (fl"5l) 
directly follows the differentiation of eq. (fl4)) utilizing the property 11401 eq. 8.2.2.30]. 

C. Proof of Proposition \J} 

When the interfering signals are from L equal power sources with i.i.d. interfering channels, 
X u {i) = X u = A/, X v (p) = X v = Xi and the characteristic coefficient X iy j (A) of A becomes |f39l , 

f j = 1,2,3.. 
X id (A) = < (24) 

[ 1 J = L 

Using the above facts the proof of the Proposition Q] immediately follows from Theorem [21 Due 
to equal power allocation, average 1st and 2nd hop SNR Ai = A 2 = A. 

D. Proof of Proposition [2] 

If the interfering channels are i.i.d. Nakagami-m faded, in Definition^ the PDF of r.v. U and 



V can be modified to, 



f v (u) = '- = e~ - , u > 1 (25) 



fv {v) = ~ L 2 m 2 e ~^> v > 1 (26) 



and 

(y - („-!) 

T (L 2 m 2 ) at 

respectively, where a u = X u /m u and ct^ = X v /m v respectively, and m is the corresponding 
Nakagami-m channel parameter. The r.v.s X and Y are exponentially distributed with hazard rate 
1/ X x and 1/ A y respectively. The CDF of F lsmR2 (7) for i.i.d. Nakagami-m interfering channels 
can be obtained by replacing X u = a u , X v = a v , j = L\m\ and q = L 2 m 2 in (fT4l) . 
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a) System model 1 b) System model 2 



Fig. 1. Interference relay network. 
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Fig. 2. Outage probability vs total SNR in system model 1 and 2. 
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Fig. 3. Outage probability vs total SNR in system model 2 with different level of INR's. 




Fig. 4. Outage probability vs total SNR in system model 2 with different number of interferers when INR=3 dB. 
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Fig. 5. Outage probability as a function of total interferers at the relay and destination when the INR of each interferer is 3 
dB and total transmit SNR varies from 20 to 60 dB. 




Interference power in dB 



Fig. 6. Outage probability as a function of interference power when the relay and the destination have 4 interferers and A to t 
is 20, 30 and 40 dB's. 
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Fig. 7. Probability density function of SINR of system model 2. 
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Fig. 8. Cumulative distribution function of SINR of system model 2. 
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Fig. 9. Outage probability vs total SNR over Nakagami-m faded interfering channels in system model 2 with different number 
of interferers. 
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Bappi Barua, Mehran Abolhasan, Daniel Franklin and Farzad Safaei 



Abstract 

Cooperative relaying improves the performance of wireless networks by forming a network of 
multiple independent virtual sources transmitting the same information as the source node. However, 
interference induced in the network reduces the performance of cooperative communications. In this 
work the statistical properties, the cumulative distribution function (CDF) and the probability density 
function (PDF) for a basic dual hop cooperative relay network with an arbitrary number of interferers 
over Rayleigh fading channels are derived. Two system models are considered: in the first system model, 
the interferers are only at the relay node; and in the second system model, interferers are both at the 
relay and the destination. This work is further extended to Nakagami-rn faded interfering channels. 
Simulation results are presented on outage probability performance to verify the theoretical analysis. 
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Signal to interference plus noise ratio, interference to noise ratio, amplify and forward, probability 
density function, cumulative distribution function. 
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I. Introduction 



Wireless channels often suffer severe performance degradation due to multipath fading and 
interference. The increase of spectral reuse in wireless systems exposes the network to ever- 
large number of interfering nodes. Ad-hoc and multiuser networks often face interference from 
other user nodes operating in multiple parallel channels. Interference can also appear in the 
network as cochannel interference. Cooperative relay communications, originally introduced 
by van der Meulen [?] has been extensively studied in [?], [?], [?], [?] without the presence 
of interference. In cooperative relay networks, cooperating nodes help the information source 
node by jointly forming a virtual antenna array and thus providing higher diversity gain and 
information reception reliability. In addition, cooperating relays can offer a very flexible extension 
of network coverage without substantial development of infrastructure. However to date, research 
efforts on the interference channels have mainly concentrated on information theoretic analysis 
[?], [?], [?], [?], [?], [?] and using assisted relay networks in [?], [?], [?], [?], [?], [?], [?], [?], 
[?], [?], [?], [?], [?], [?], [?]. Statistical modeling for single hop multiple antenna systems with 
active interferers has been performed for maximal ratio combining (MRC) receivers in [?], [?], 
[?], [?], [?] and using optimum combining receivers in [?], [?], [?], [?]. 

Interference in cooperative relay networks has been extensively studied in previous publications 
[?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?]• Authors of [?] and [?] consider a scenario 
where only the relay node experiences interference in interference limited networks and the 
receiver node remains interference free. The relays here can estimate the instantaneous CSI of 
interfering channels to scale the gain. A fixed gain AF relay with interference limited destination 
is considered in [?]. However, the assumption that the AF relay gain parameter includes the 
instantaneous or average channel information of interfering channels (as assumed in [?], [?], 
[?], [?], [?], [?], [?]) requires additional computational capability at the relaying node, and in 
certain cases where the interfering signals are not known to the relay a priori, the technique 
can not be applied. Furthermore, in many previous results authors consider interference limited 
cooperative networks [?], [?], [?], [?], [?], [?]. Outage performance of a dual hop network has 
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been studied using a fixed gain relay in [?] and hypothetical gain AF relay in [?] with an 
arbitrary number of interferers. The system is assumed to be an interference dominated network 
where noise power is negligible compared to interfering signal power. Closed form PDF and 
CDF of signal to interference ration (SIR) are derived. Avoiding the noise in receiver systems 
usually provides an opportunity to express some of the integrals used to obtain the CDF and 
the PDF in very standard formats. Outage probability using DF relays in Nakagami-m fading 
channels is studied in [?]. The authors consider that the destination faces a negligible amount 
of interference, but non-negligible noise. The PDF and CDF of corresponding receiver SNR is 
derived. In [?], the authors have studied transmission techniques designed to manage interference 
using an information-theoretic approach. Bit error rate (BER) performance has been presented in 
[?] for a multi-user network in which multiple sources communicate with one receiver node via 
a common relay in a dual hop network. A different approach has also been proposed to cancel 
the interference at the relay node. In [?], the authors have investigated the effect of cooperation 
in an interference-limited system utilizing a decode-and-forward (DF) relaying protocol over 
Rayleigh fading channels. It divides the cooperating nodes into different cooperating regions 
and quantifies the relation between cooperative region radius and the interference level. It also 
analyzes the network sum rate optimization problem for a cooperative region radius. 

In this work, we derive the closed form expressions of statistical properties, the CDF and 
the PDF of signal to noise plus interference ratio (SINR) of an interference relay network that 
is operating under the influence of arbitrary number of interferes and noise. An amplify-and- 
forward (AF) relaying protocol is adopted due to its simplicity and ease of deployment at relay 
nodes. Lastly numerical results on the outage probability performance are presented to verify 
the theoretical analysis via Monte-Carlo simulations. 

II. System and Channel Model 

Consider a single source-destination pair communicating via a single antenna relay without any 
direct link. We will denote source-relay and relay-destination links as S-R and R-D respectively. 
A half duplex AF protocol has been considered over independent and non identically distributed 
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(i.n.i.d.) Rayleigh fading channels. In this paper, two different system models are investigated: 
system model 1 (SM 1), in which the interferers are only at the relay node; and system model 
2 (SM 2), where interferers affect both the relay and destination nodes. In both models, all the 
interfering channels are i.n.i.d. Rayleigh faded. The destination is assumed to have full channel 
state information (CSI) of the two main channels, S-R and R-D, while the relay has full CSI of 
the S-R channel only. The source and relay have no CSI of forwarding transmitting channels. 
None of the nodes, source (S), relay (R) and destination (D) possess information about the 
interfering channels. 

The instantaneous and average signal to noise ratio (SNR) of 1st and 2nd hops are denoted as 
7i — P][h 'i and A, = respectively, where j 6 {1,2}, Pis the corresponding source and relay 

i i 

power; hi and f2j are the instantaneous and average channel gain of the ith hop respectively^ and 
of is the one sided additive white Gaussian noise (AWGN) power at relay or destination node, 
i.e. i E {R, D}. We assume the total power of the network (source and relay) is constrained to 
P tot . This total power is split between the source and the relay by a power sharing coefficient 
C £ (0, 1] such that the source and the relay powers are given by, Ps = (1— C)Pot and Pr = (Pot 
respectively. Thus if equal power sharing protocol is adopted, £ will be 0.5. 

Let there be a total of L interferers in the system, and define an interferer set J, the set of all 

interfering source nodes. For example, any interferer // £ 3, where I E £, £ = {1,2,3 , L}. 

The elements of interfering channel row vectors hf j E C L represent the corresponding interfer- 
ence channels from the source element nodes of the interfering signal source vectors xj $ E C L , 
i E {R, D}a Furthermore, we assume the interferer Ii has individual transmit power Pj t . Thus 
the instantaneous and the average interference power to noise ratio (INR) for any interferer // 

P \h | 2 

is 7i,z — — / and \j i = Pj, '^ 7,i respectively, hij and of are the fading channel gain from 
the interfering source to the node i and the noise power at node i respectively, / E C and 
i E {R, D}. ttjj is the average interfering channel gain, Qj i = E [|/i^| 2 ] . Throughout this paper 



'The average channel gain, Q,i is in fact the statistical average of the squared instantaneous channel gain hi, i.e. 0; = E [|/ii| 2 ] . 
2 C L denotes a L-dimensional complex vector. 
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the term INR will be used to indicate individual interferer signal power to noise ratio unless 
otherwise specified. 



A. System Model 1 : Interference at Relay 

In the network model as shown in Fig. [H interference occurs only at the relay node. The 
received signal at the relay node R is 

y R = hix s + h^ixj^ + rii (1) 

where, hj x G C L and xj .1 G C L are the row vectors for the interference channels and the 
corresponding interfering signal sources respectively. The notation (.) T represents the transpose 
operation on a vector. The signal at the destination node is 

y D = Gh 2 h 1 x s + Gh 2 h Ij iy:J 1 + Gh 2 n x + n 2 (2) 

where G is the AF relay gain, ni ~ CM (0, of) and n 2 ~ CM (0, a 2 ) are additive white Gaussian 
noise at the relay and destination respectively^ At the receiver node, the signal to interference 
plus noise ratio (SINR) will be 

G 2 \hi\ 2 \h 2 \ 2 P s 

7SINRl " G^h^j^ih^ + G^faj + al ( } 

where the diagonal matrix Ej = E |x| iX/^j is composed of interference signal powers. With 
hypothetical AF relay gain G 2 = p ^^ a , the end-to-end SINR 



is [?fl 



7i72 , A N 

TSINRi = : : (4) 

7i + 72 + 727/ 

where 71 and 7 2 are the instantaneous SNRs for S-R and R-D links respectively. 7/ is the total 
interference to noise ratio at the relay node, which is literally the sum of all individual INRs. 

Note that due to the Rayleigh fading assumption, the first and second hop SNR ji and 72 are 
exponentially distributed with mean parameter Ai and A2 respectively. Later, it will be shown that 

3 CA/" (A, cr 2 ) denotes a circularly symmetric complex Gaussian random variable with mean A and variance a 2 . 
4 A CSI assisted relay gain G 2 = ^ +(y >± proposed in [?] render the end-to-end SINR as 7sinr = 71+7 ^_^ 7f +1 ■ Eq- © 
proposes a tight upper bound on the CSI assisted SINR 7sinr. 
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under an equal power sharing protocol at the source and relay nodes, the system performance 
metric of system model 1 does not vary if the interferers are switched from the relay to the 
destination node. 



B. System Model 2: Interference at Relay and Destination 

In this network model, interference occurs both at the relay and at the destination nodes. 
We let there are L x interferers at the relay node and L 2 interferers at the destination. All the 
interfering channels are independent and non-identically distributed. In this case, the received 
signal at the destination is 

y D = Gh 2 h x x s + Gh 2 h I:1 xJ 1 + h /j2 x^ 2 + Gh 2 n x + n 2 (5) 

where x/ i e C Ll and x/ 2 £ C L2 are the vectors with interference sources for the relay and 
receiver nodes respectively, and h/i e C Ll and h/ 2 € C Li are the corresponding fading channels 
from interferers to the relay and receiver nodes respectively. Again, as for SM 1, G is the AF 
relay gain, n x ~ CM (0, a\) and n 2 ~ CM (0, of) are AWGN at the relay and the destination 
respectively. Thus the SINR with arbitrary relay gain is given by 

G^N'h^E^h^ + h /i2 S /)2 h} i2 + G*\h 2 \ z af + a\ 

S/ ! = E |x| iX/ ij and X/ 2 = E |x| 2 X/ 2 | are diagonal matrices of the transmission powers 
of interfering signals at the relay node and the destination respectively. Applying the similar 
hypothetical relay gain as used in SM 1 at the AF relay, the receiver SINR is given by 

7i72 ( ~. 

7SINR 2 = : : : (/) 

7i + 72 + 7i7/,2 + 727/,i 

where 7/1 and 7/ j2 are the total interference to noise ratio at the relay node and the destination 

respectively. 



7SINR 2 - _ ol , l2 , -f — ^ 71 . ,2 o . o 



III. Statistical Analysis 
In the following sections statistical theorems related to the network of interest are developed. 
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Definition 1: Let random variables X, Y and U be statistically independent, where X 
£ {j^j an d Y ~ £ {j^j J3 Random variable U — J2 Ui + 1, where Ui ~ £ (x - ) are i-n.i.d. expo 



C(A)r»(A) 



nentially distributed random variables and U has the PDF, fy {u) = Yl Yl Xl,J ^^ u j 1 - > — e 

i=i j=i r w A «<i> 

u > 1 [?]. A r.v. W is defined as 



n 



W = XY (8) 

Theorem 1 (CDF and PDF): The cumulative distribution function (CDF) and probability den- 
sity function (PDF) of random variable W are given respectively as 



C(A) Ti(A) 

F w («,) = 1 - £ jfy (A)e" 

i=l j = l 



A, 



oo k+1 

*+EE 



r (j) (X x Xy) + k\(k + 1)! VA + w\ u 



X x + wX u (i) 

2T (j +n)Un 



k=0 n=0 



k + 1 

n 



w 



\/X X Xy 



+ G 



1,3 
.3.2 



/ X X X U ^ 


l-J-n,l,l\ 


\X X + w\ u (i) 


1,0 / 



(9) 



C(A)7i(A) 

fw h = E E x « ( A ) e 
i=i j=i 

OO fc + l ^ 

-EE 



An 



As + wX u (i) 



k=0 n=0 



w 



2fc+l 



T(j)(X x X y ) k+1 k\(k + l)\ \X x + wX u{l) 



X x X u (i) 



_A X Ay A x + wX u (i) 

2T (j + n) + Zi (w, A) <( 2A; + 2 



tw (A* + Ay) (j + n) wA„(i) I wX u{i) ,, : . 

+ ~ — : : ^4,3 



X x X y 



X x + wX u {i) X x + wX u (i) 



X x X u (i) 
X x + wX u (i) 



1-j-n, 1,1,0 
1,1,0 



(10) 



5 £ (x - ) denotes an exponential distribution with a hazard rate j-. 

6 where A = diag (A U1 , \ U2 , A„ L ). £ (A) is the number of distinct diagonal elements of A and X u n\ > A u ( 2 > > > 

Au(c(A)) are the distinct diagonal elements in decreasing order, n (A) is the multiplicity of \ u <i) and Xi t j (A) is the (i, j')th 
characteristic coefficient of A [?]. 
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where ^ = \ V (* + 1) + \ 1> (k + 2), and X l (w,X) ± G\% 



1 - j - n, 1, 1 
1,0 



+ 



Up 



and ijj (x) are the Meijer-G function and the 



Euler psi function defined as [?, eq. 8.2.1.1] and [?, eq. 8.360.1] respectively. 

Note that, exchanging X and Y in Definition 1, represents the switching of the interference 
from relay to destination node in ©. Hence, the CDF and PDF equations of SMI are also 
applicable to the general case when either the relay or destination is subject to interference. In 
addition, the performance metric will be the same if an equal power sharing is adopted between 
the source and relay nodes, that is, when ( = 0.5. 

Proof: See Appendix [A] ■ 

Definition 2: Let random variables X, Y, U and V are statistically independent where X ~ 
£ (i) and Y ~ £ (i) . Random variables U = £ U t + 1 and V = £ V x + 1, U t and Vj are 
i.n.i.d. exponentially distributed random variables with PDFs [?] 



fuM= E E F7T^( u - i r le " A " <4> . 

i=l j=l 1 UJ 



C(A 2 )r,(A 2 ) 

/v W=EE WTZT-b - v > 1 (12) 



respectively. We define a r.v. Z such that 



Z = XY (13) 
Xf/ + YV 



Theorem 2 (CDF and PDF): The cumulative distribution function (CDF) and probability den- 
sity function (PDF) of random variable Z are given respectively as 
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C(A 1 )r 8 (A 1 )C(A 2 )r p (A 2 ) 

^W = 1 "EEEE X m (Ai) X M (A 2 ) e 
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AyA u (;) 



fc=0 ra=0 r=0 
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(14) 
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1 — g — n, 1, 1, 
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H : : 



Ay + zX u (i) 4 ' 3 ^ Ay + ^A u (j) 

where, 

X 2 (z, A) = 2T (g + n) T (j + r) <j In 



1 -j -r, 1,1,0 



1,1,0 



(15) 



+ T(g + n) Gg' 2 



1,3 I 



Proof: See Appendix [B 



1 — j — r, 1, 1 
1,0 



1 — g — n, 1, 1 
1,0 
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Proposition 1: Consider a system with L interferers in both relay and destination nodes where 
the interfering channels are i.i.d. Rayleigh faded. Adopting an equal power allocation protocol 
for source and relay node the cumulative distribution function of 7sinr 2 in © can be written as 



F (V) = 1 - e" 27/A ( - 



2L 



T(L)' 



-2 7 /a 



00 fe+i fe+i 

k=0 n=0 r=0 



k + l\ (k + 1 



n \ r 



n+r 2k+2 



X 



\ n+r 

Aj 7 



A 



2L+n+r 



2r (L + n)T(L + r) {in^-^} 



+r [L + r) 



AA, 



A + 7A7 



1-L-n, 1,1 



1,0 



+ r(i + n)G ,3|_^_ 



1 -L-r, 1,1 



1,0 



(16) 



Proof: See Appendix | 



IV. Interference at Relay: I.I.D. NAKAGAMi-m Interferers 

Proposition 2: Suppose the interfering channels are i.i.d. Nakagami-m distributed while the 
main channels (S-R and R-D) are Rayleigh faded. In this case the CDF of 7sinr 2 can be written 

as 



7SINR 2 



( 7 ) = l-e" 
00 fc+i fc+i 
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Aiaj,2 
Ai + 7011,2 
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n r ;r(m 1 L 1 )r(m 2 L 2 )A;!(A; + l)! VA 2 + 7a 7 ,i 



m 2 L 2 



2T [m\L\ + r) T (m 2 L 2 + n) < In 



7 
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A 2 aj,i 
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1,0 

1 — rriiLi — r, 1, 1 
1,0 



(17) 
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where ai t i = Xi^/rrii, i E {1, 2}. A/,1 and A/^ are the average INR at the relay and the destination 
respectively; similarly, mi and rri2 are i.i.d. Nakagami-m parameter at the relay and destination 
respectively. The average SNR of the 1st and 2nd hop are Ai and A 2 respectively. 

Proof: See Appendix iDl ■ 

V. Outage Probability 

In this section the derived analytical results are used to the investigate wireless network 
outage probability performance. We define outage probability as the probability that the in- 
stantaneous receiver SINR falls below a predefined threshold value of SINR 7th. We consider 
7th = P (2 MR — l), where p varies from 1 to 6.4 depending on the degree of coding, M be 
the number of hops and R the data rate in bits/s/Hz [?]. The CDF equations ®, CLU) and (fTTT) 
may be used to evaluate the outage probabilities in SM 1 and SM 2 with an arbitrary number 
of interferers and interfering powers. 

VI. Numerical Analysis 

In this section, numerical results on the derived equations are presented and compared with 
Monte-Carlo simulations. Due to lack of transmitter CSI we assume the source and the relay 
evenly share the total system power P tot . Furthermore, for fair comparison with relevant studies 
the average channel gain of all hops in this network are set to unity and the noise variance at 
the relay and destination nodes are set to a 2 . We consider 7 th = 3 which achieve a data rate 
R — 1 in a dual hop network unless otherwise specified. The maximum value of k = 100 for 
sum terms in © and (fl4l) is used, which is sufficient to match the analytical result with the 
simulation data. The number of interferers at the relay and destination are assumed to be same 
(Li = L 2 ) in SM 2 networks, and all the interferers have equal transmit power. Throughout this 
section we will use total average transmit SNR, A to t — -Ptot/c 2 as a function in performance of 
SINR plots. 

Fig. [2] shows the outage probability (OP) as a function of total transmit SNR for a total of 
{4, 8, 16} interferers in SM 1 and SM 2 with INR fixed at 3 dB. All wireless channels (main 
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channels and the interfering channels) are assumed to be Rayleigh faded. As expected, outage 
probability increases with an increase in the number of interferers. However, when SM 1 has 
the same number of interferers as the total number of interferers for SM 2, SMI shows slightly 
better performance compared to SM2 at a moderately low SINR region only. To study the OP 
performance with different interfering powers Fig. [3] shows outage probability as a function of 
A to t for a system with 4 interferers at the relay and 4 interferers at the destination. The INR 
values {3, 6, 9, 12} dB are considered. Interestingly, the figure suggests that every 3 dB increase 
in average INR values requires a subsequent increase of 3 dB total SNR on average to ensure 
an OP of around 1CT 3 . In contrast, Fig. @] shows OP in a fixed 3 dB INR system in SM 2 with 
a total of {4, 8, 12, 16} number of interferers at the relay and destination respectively. To see 
how the interference system behaves if one of the nodes experiences most of the interference, 
Fig. H] shows OP with a large L 2 and L\ = 4. For very high interference levels e.g. L 2 = 10000 
interferers at the destination node, a A tot = 77 dB is required to achieve an outage probability 
of 10~ 3 compared to 50dB when L 2 = 4. 

Fig. [5] shows outage probability with number of interferers L in relay and destinations where 
Li = L 2 = L and each interferer has an INR of 3dB. The figure shows how the outage 
probability increases exponentially as the number of interferers increases. The plot includes 
graphs for transmit SNRs of 20, 30, 40, 50 and 60 dB. A very interesting phenomenon observed 
is that at 60dB of total transmit SNR the system can support up to 40 interferers at the relay 
and destination, ensuring an outage probability of 10~ 3 , however, with a reduction of only lOdB 
in total transmit SNR the performance reduces drastically and it can support only 4 interferers 
at the relay and destination for the same outage probability. 

Fig. [6] shows the outage probability as a function of interference power when the relay 
and destination are both subject to 4 interferers. The figure compares the outage probability 
of an interference limited system (dashed line) studied in [?] to a system where both noise 
and interference are present. To express total transmit power and interference power in dB we 
normalize the terms considering noise power equal to 1 unit. Total transmission power of 20, 



BARUA et al : OUTAGE PROBABILITY ANALYSIS OF DUAL HOP RELAY NETWORKS IN PRESENCE OF INTERFERENCE 13 

30 and 40 dB are considered in the figure and the power of each interferer is rated from -30dB 
to +30dB. Figure shows that the outage probability the noise plus interference network will be 
similar to an interference limited network when the transmission power of each interferer is more 
than 10 dB. However, performance will vary if the number of interferers in the network changes. 
The figure also illustrates how noise dominates outage probability metric when interference power 
is low. 

The PDF of 7sinr 2 is presented in Fig. [7] for two different A tot , 20 dB and 25 dB, when the 
individual interferer INRs are {3, 6, 9} dB. The figure implies that a lower INR per interferer 
increases the probability of higher output SINR at the receiver and similarly when the INR is 
constant, an increase in average total signal power results in a higher output SINR. In addition, 
Fig. [8] shows the CDF of 7sinr 2 as a function of 7 th in dB for two A tot groups, 20 dB and 30 
dB. Each group contains plots for 3 INRs 3, 6 and 9 dB. It reveals that higher average total 
SNR and lower INR per interferer decreases the probability of outage. 

Finally, Fig. [9] presents a plot over Nakagami-m faded interfering channels. The outage 
probability is plotted as a function of A to t for different number of interferers and different 
Nakagami channel parameter m in SM 2. Surprisingly, varying the Nakagami m parameter does 
not result in any significant effects on outage performance if the number of interferers and INR 
remain constant in the network. 

VII. Conclusion 

In this work, we derive the cumulative distribution function and the probability density function 
of a dual hop interference relay network with arbitrary number of interferers, where the main 
channels and the interfering channels are i.n.i.d. Rayleigh faded. Later the analysis is extended 
for i.i.d. Nakagami-m faded interfering channels. This analysis has been performed for two 
system models; in SM 1 interference is only at the relay node while in SM 2, interference 
is at both the relay and destination nodes. The derived CDF and PDF expressions for SM 1 
can also be used for analyzing an interference network where interference occurs only at the 
destination node. Numerical results for outage probability performance for different network 
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parameter configurations are presented. The Monte-Carlo simulations show an exact match with 
the analytical expressions. Results indicate that in a total INR constrained interference network, 
the number of interferers do not affect the performance of the system in a large scale. Similarly, 
if the interfering channels are Nakagami-m distributed, changing in Nakagami-parameter m does 
not affect the performance of the system if the number of interferers and the INR per interferer 
remain constant. 



Appendix 

A. Proof of Theorem [7J 

According to the PDF of X, Y and U as defined in Definition 1, the random variables X, Y 
and U are nonnegative, thus F w (w) = for w < 0. For w>0we have 

r XY 

= P {X (Y - w) < wYU} 

wyu 

poo poo r~_ — 

= 1 - e~ w/x « + / / I fx 0) h (V) fu (u) dxdydu (18) 
Ji i» Jo 



The two inner integrals in (1181) can easily be solved by using [?, eq. 3.351.3] and [?, 
eq. 3.471.9] and after some manipulation the CDF of W is obtained, 



C(A)r l (A) , . p -w/\ v roo 

i= i j= i 1 UJ A u(i)V A * A v 



e '™" x ' Ki ( 2w \[S) du (19) 



where K v (x) is the z/th order modified Bessel function of the second kind. Using [?, eq. 8.446] 
we expand the first order modified Bessel function of second kind, then [?, eq. 8.4.6.5], [?, 
eq. 3.40.1.1] and [?, eq. 3.351.3] are applied to solve the related integrals. After some manipu- 
lation we arrive at the desired result ©. The PDF of W follows directly from the differentiation 
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of the CDF of W in eq. © w.r.t. w by using [?, eq. 8.2.2.30], 



d_ 

dx 



30 q I 



tj p+l,q+\ 
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dp, -a \ 








l-a,b q J 



(20) 



B. Proof of Theorem \2\ 

Following a similar argument as used in the proof of Theorem \T\ the random variables X, Y, 
U and V are nonnegative, and thus F z (z) = for z < 0. For z > 0, 



F z (z) = P 



XY 



< z 



XU + YV 
¥{X(Y -zU) < zYV} 
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(21) 



fx (x) f Y (y) fv {v) fu (u) dxdydvdu 



A T 
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V 

- T 

— J-uvyx 



The first part of this expression, integral I uy can be solved by using [?, eq. 3.351.3], resulting 
in 



(22) 



C(Ai) ri(Ai) , sj 

In I uvyx , the first two integrals for variables x and y can be solved using a similar approach, via 
[?, eq. 3.351.3], [?, eq. 3.471.9], [?, eq. 8.446], [?, eq. 8.4.6.5] and [?, eq. 3.40.1.1]. Substituting 



(1221) in (|2lT) and after some manipulation we have the CDF of Z 
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(23) 
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where ip k = \ (k + 1) + \ if) (k + 2). Now using [?, eq. 3.351.3], [?, eq. 8.4.6.5] and [?, 
eq. 3.40.1.1] it is possible to solve (l23l) . which then results in the desired eq. (fl4l) . Eq (fl"5T) 
directly follows the differentiation of eq. (fi4")) utilizing the property [?, eq. 8.2.2.30]. 

C. Proof of Proposition \J} 

When the interfering signals are from L equal power sources with i.i.d. interfering channels, 
X u {i) = X u = A/, X v (p) = X v = A/ and the characteristic coefficient X it j (A) of A becomes [?], 

f j = 1,2,3.. 
X id (A) = < (24) 

[ 1 3 = L 

Using the above facts the proof of the Proposition Q] immediately follows from Theorem [2l Due 
to equal power allocation, average 1st and 2nd hop SNR A x = A 2 = A. 

D. Proof of Proposition [2] 

If the interfering channels are i.i.d. Nakagami-m faded, in Definition^ the PDF of r.v. U and 
V can be modified to, 

f v ( u ) = A 1 e" - , u > 1 (25) 

and 

fv «) = ^ \ W ~ ^ > 1 (26) 

F (L 2 m2)a v 2 

respectively, where a u = X u /m u and a v = X v /m v respectively, and m is the corresponding 
Nakagami-m channel parameter. The r.v.s X and Y are exponentially distributed with hazard rate 
1/ X x and l/Xy respectively. The CDF of -^ysinr (7) f° r Nakagami-m interfering channels 
can be obtained by replacing X u = a u , X v = a v , j = L x m\ and q = L 2 m 2 in (fT4l) . 
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a) System model 1 b) System model 2 



Fig. 1. Interference relay network. 




20 30 40 50 



Kx ( dB ) 

Fig. 2. Outage probability vs total SNR in system model 1 and 2. 
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20 30 40 50 60 



Fig. 3. Outage probability vs total SNR in system model 2 with different level of INR's. 




Fig. 4. Outage probability vs total SNR in system model 2 with different number of interferers when INR=3 dB. 
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Fig. 5. Outage probability as a function of total interferers at the relay and destination when the INR of each interferer is 3 
dB and total transmit SNR varies from 20 to 60 dB. 




Interference power in dB 



Fig. 6. Outage probability as a function of interference power when the relay and the destination have 4 interferers and A to t 
is 20, 30 and 40 dB's. 
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Fig. 7. Probability density function of SINR of system model 2. 




SINR Threshold (dB) 



Fig. 8. Cumulative distribution function of SINR of system model 2. 
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Fig. 9. Outage probability vs total SNR over Nakagami-m faded interfering channels in system model 2 with different number 
of interferers. 



